2009 Fermi Symposium, Washington, D.C., Nov. 2-5 


1 


4U 1909+07: a well-hidden pearl 
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We present the first detailed spectral and timing analysis of the High Mass X-ray Binary (HMXB) 
4U 1909+07 with INTEGRAL and RXTE. 4U 1909+07 is detected with an average of 2.4 cps in 
ISGRI, but shows flares up to ~50 cps. The system shows a pulse period of 605 s, but we found that 
the period changes erratically around this value. The pulse profile is extremely energy dependent: 
while it shows a double peaked structure at low energies, the secondary pulse decreases rapidly with 
increasing energy and above 20keV only the primary pulse is visible. This evolution is consistent 
between PCA, HEXTE and ISGRI. We find that the phase averaged spectrum can be well fitted 
with a photoabsorbed power law with a cutoff at. high energies and a blackbody component. To 
investigate the peculiar pulse profile, we performed phase resolved spectral analysis. We find that a 
change in the cutoff energy is required to fit the changing spectrum of the different pulse phases. 
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I. INTRODUCTION 

The X-ray source 4U 1909+07 (or XI 908+07) was 
detected with the Uhuru satellite [1] as 3U 1912+07 
and filed in the fourth Uhuru catalog as 4U 1909+07 
[2], It has been seen with most X-ray instruments ever 
since. The exact nature of the system was. however, 
not clear for almost 30 years. In 2000, Wen et al. [3] 
analyzed RXTE ASM data and found a stable 4.4 days 
period, which was interpreted as the orbital period of a 
binary orbit. Due to the high absorption, however, no 
optical counterpart could be identified. Using RXTE- 
PCA data, Levine et al. [4] found a second, shorter 
period of ~605s in the X-ray flux, explained as the 
pulse period of a slowly rotating neutron star. Using 
Doppler delays curve Levine et al. [4] could refine the 
binary orbit parameters. They estimated the mass 
of the companion star to be M* = 9-3 1 M 0 using 
a canonical mass of 1.4 Mg for the neutron star and 
the radius of the companion star to be /?* < 22 Rq. 
One year later Morel and Grosdidier [5] detected an 
OB star in the near infrared at the location of the 
X-ray source, thus confirming that the system is a. 
High Mass X-ray Binary (HMXB). The distance of 
the system was estimated to be 7kpc [5]. Prior to 
this discovery, Levine et al. [4] argued that the com- 
panion star could be a Wolf-Rayet star, which would 
make this system a possible progenitor system to a 
neutron star-black hole binary. This possibility, al- 
though intriguing, can be ruled out after the discov- 
ery of Morel and Grosdidier [5]. Although the system 
shows no eclipse the X-ray flux is- still strongly orbital 
dependent, as seen in ASM data folded onto the or- 
bital period [4]. Levine et al. [4] also analyzed orbital 
resolved PCA spectra and found that around orbital 


phase 1 the photoabsorption increases dramatically by 
a factor of at least 3, explaining the decreased ASM 
flux. This increase in absorption can be very well de- 
scribed by a spherical wind model and an inclination 
of 54° <i< 70°, depending on the parameters of the 
wind model. Like most wind-fed HMXB with a neu- 
tron star, 4U 1909+07 shows a constant , irregular flar- 
ing behavior, superimposed on the regular pulsations. 
These flares are most probably due to inhomogeneities 
in the wind of the donor star, which lead to different 
accretion rates. In these proceedings we present the 
data and reduction methods in Sect. II. In Sect. Ill 
we analyze the pulse period evolution and the pulse 
profiles and in Sect. IV we perform phase averaged 
and phase resolved spectroscopy. We summarize and 
discuss our results in Sect. V. 


II. DATA AND OBSERVATIONS 

For our study of the 4U 1909+07 system we used 
all available public data of the X-ray satellites IN- 
TEGRAL [6] and RXTE, INTEGRAL features the 
detector ISGRI. which is a coded mask instrument 
sensitive in the 15keV - lOMeV energy range and 
part of the Imager on Board the Integral Satellite 
[IBIS: 7, 81. Even though no pointed observations on 
4U 1909+07 were performed with INTEGRAL, thanks 
to ISGRI’s large field of view of almost 30° x 30° 
more than 6 Msec of off-axis data exist. A huge part 
stems from the core program on the galactic center 
performed in the early AOs [9]. Data from later AOs 
was taken from the extended monitoring campaign of 
the mieroquasar GRS 1915+105 [10]. Twelve sources 
with a detection significance > 6 a were found in an 
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ISGRI mosaic consisting of the images of all ScWs, 
with 4U 1909+07 being the fourth brightest source. 
The data were extracted using the standard pipeline 
of the Offline Science Analysis (OSA) 7.0 for spectra 
and images and using iiJLight to obtain lightcurves 
with a higher temporal resolution. We extracted the 
lightcurve with a resolution of 20 sec as a trade-off 
between good signal-to-nose ratio and good time res- 
olution to measure the pulse period. 

For phase resolved spectroscopy we used RXTE ob- 
servations carried out in early 2003 which were part 
of the analysis done by Levine et al. [4], These data 
had a total exposure of 196ksec. Both major instru- 
ments of RXTE were fully operational at that time, 
the Proportional Counter Array [PC A; 11], sensitive 
in the 2-60keV energy band and the High Energy X- 
Ray Timing Experiment [HEXTE; 12], sensitive be- 
tween 1 5-250 keV. Both instruments have a large ef- 
fective area and a high temporal resolution, but no 
imaging capabilities. RXTE data were reduced using 
the standard HEASARC software (v. 6.6), obtaining 
lightcurves of the PCA with 1 sec resolution and pulse 
phase resolved spectra for PCA and HEXTE in the 
4.5-150keV range. 


III. TIMING ANALYSIS 
A. Pulse period evolution 

As the pulse period of the system was only discov- 
ered in 2004 using RXTE data from 2000 [4] , no his- 
toric data on the pulse period evolution exist. Be- 
tween 2000 and 2004, the pulse stayed at a value of 
P « 605.68s [4], The archival INTEGRAL data now 
provide an optimal basis to perform many more mea- 
surements of the pulse period between 2003 and 2008 
and to follow its evolution. For this analysis we split- 
ted the lightcurve into segments between 300 ksec and 
SOOksec length, providing a good balance between ac- 
curate pulse period determination and temporal res- 
olution, and used epoch folding [13] to determine the 
pulse period for each of these segments individually. 
4U 1909+07 shows strong pulse-to-pulse variations, in 
fact, the pulse period is not visible in the ISGRI 
lightcurves to the naked eye, see, e.g., Fig. 1. Only 
after folding 500-1000 pulse periods, a stable pulse 
profile arises which allows for a reliable pulse period 
determination. The emerging profile is, however, re- 
markably stable over our data set. A similar effect is 
seen in Vela X-l, although the pulsed fraction and rhe 
overall luminosity are distinctly higher in that source 
[14]. 

The pulse period between 2001 and 2007 is shown 
in Fig. 2. Albeit two trends could be picked out, a 
spin-up trend between March 2003 and April 2006 and 
a spin-down trend from April 2006 to October 2007 
the overall behavior is consistent with a random walk 
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FIG. 1. Small part of the lightcurve in the 20-40 keV en- 
ergy band as measured with ISGRI. The red line shows 
the folded pulse profile over the whole data set. 

like behavior. To confirm this observation, we imple- 
mented the algorithm proposed by de Kool and Anzer 
[15]. This algorithm evaluates the relative change in 
pulse period over different time intervals 5t, between 
single measurements. For a perfect random walk, the 
result will be a straight line with a slope of 0.5 in 
logdu; — logdt space, where ui is the angular velocity: 
u) = 27 t/P. As only very few measurements of the 
pulse period of 4U 1909+07 were available, the out- 
come of the calculation is relatively coarse. Nonethe- 
less, it is clear from Fig. 3 that the data follow a 
straight line. Superimposed on the plot is a line with 
slope of 0.5, shifted in y-direction to fit the data. The 
uncertainties to the data are not only due to uncer- 
tainties in the estimate of the pulse period, but are 
also due to the uneven and coarse sampling. To esti- 
mate this effect, we choose a Monte Carlo approach 
and simulated many mock-up pulse period evolutions 
which followed a perfect random walk and sampled 
these with the same rate as the real data [see 15]. 
The standard deviation in each bin from these simu- 
lations gives an estimate of the uncertainty in each bin 
and is roughly of the same order of magnitude as the 
uncertainties from the pulse period determination. 

B. Pulsprofiles 

As shown by Levine et al. [4], the 3.7 - 17keV pulse 
profile shows two distinct peaks, with the second one 
being slightly broader than the first one and having a 
complex shape with two subpeaks at phase 4> = 0.85 
and o = 1.1. But the pulse profile of 4U 1909+07 
shows a remarkable and strong energy dependence. 
To follow this energy dependence of the pulse profile, 
we extracted pulse profiles from PCA data in differ- 
ent energy bands, using a period of P = 604.685 s and 
32 phase bins. Thanks to the large effective area of 
the instrument, very good pulse profiles could be ex- 
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FIG. 2. Evolution of the pulse period over the last years. 
The historic RXTE data points of Levine et al. [4] are 
shown as red crosses, the new measurements obtained with 
INTEGRAL are shown as stars. 
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FIG. 4. Energy resolved pulse profiles with RXTE PCA 
( a)~c )) and INTEGRAL ISGRI (d)). The profiles are 
shown twice for clarity. Note that the RXTE and IN- 
TEGRAL profiles are not phase aligned. 

of the data. The trend seen in RXTE pulse profiles 
is continued, and the secondary peak has completely 
vanished, where as the primary is very clearly seen 
(Fig. 4d). 
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FIG. 3. Plot of the pulse period evolution in log Sui — log 5t 
space, as proposed by de Kool and Anzer [15]. A ideal 
random walk process u'ould show a slope of 0.5 in this 
plot. The data follow this slope closely, as shown by the 
red straight line wfith a slope of 0.5. 


tracted in 30 narrow energy bands. Three examples 
are shown in the upper panels of Fig. 4. The pulse 
profiles show a smooth transition from a two peaked 
profile at low energies to a single peak profile at high 
energies. At very low frequencies the secondary peak 
is broader and stronger than the first one (Fig. 4a). 
This behavior levels out with increasing energy and 
the two peaks become more separated (Fig. 4b). The 
relative power of the secondary peak keeps declining 
with increasing energy and is very small compared to 
the first one at energies around 20keV (Fig. 4c). The 
deep minimum around phase 0.3 is not energy depen- 
dent. For the very high energies we extracted the pulse 
profile from INTEGRAL data of the average data of 
the first 100 days of measurement in 2003. We used 
the same epoch as for the PCA analysis, but a pe- 
riod of P — 604.747 s, as determined in our analysis 


IV. SPECTRAL ANALYSIS 
A. Phase averaged spectrum 

The pulse phased averaged spectrum of 4U 1909+07 
show r s a shape similar to many other neutron star 
sources, consisting of a powerlaw continuum atten- 
uated by photoabsorption at low- energies and by a 
turnover at high energies. Levine et al. [4] mod- 
elled the spectrum using a cutoff-powerlaw and a 
bremsstrahlung model and found that both fit the 
data equally well. The turnover at high energies 
is modeled phenomenologically with different models 
like cutoffpl or highecut, where in the latter the 
cutoff is only applied after a certain energy E cut . More 
sophisticated ones include fdcut [16] and npex [17]. 
For our spectral analysis we discarded all observa- 
tions between orbital phase 0.88 < <j> or b < 0.12, where 
the Ah is dramatically increased [4]. To this dataset 
we fitted the above mentioned models to RXTE and 
INTEGRAL data simultaneously using XSPEC 11.3. 
We found that all models give equally good fits re- 


eConf C091122 





2009 Fermi Symposium., Washington, D.C., Nov. 2-5 



FIG. 5. a) Combined RXTE PCA (red) and HEXTE 
(blue) and INTEGRAL ISGRI (green) spectra. The best 
fit cutoffpl-model is also shown, b) Residuals of the best 
fit model. 

garding the reduced x 2 value. In the fdcut model, 
however, the cutoff energy is fitted to values < 1 ke V, 
thus effectively reducing the model to a cutoffpl. 
Some of the most important spectral parameters of 
all other models are shown in Tab. I. Independent 
of the applied continuum models, a soft excess below 
lOkeV is evident, which can be very well modelled us- 
ing a blackbody with a temperature of kT a:1.4keV. 
Additionally strong evidence for an iron line close to 
the energy of neutral iron of 6.4 keV is found. The iron 
line was modelled with a Gaussian with an equivalent 
width of ~250eV. Fig. 5 shows the spectrum and the 
best fit cutoffpl with its residuals. 


B. Phase resolved spectra 

We divided the PCA and HEXTE data into 7 phase- 
bins. as indicated in Fig. 6a). The INTEGRAL ISGRI 
data did not provide high enough statistics for high 
resolution phase resolved spectroscopy. To all phase 
bins we fitted a cutoffpl with black body and a npex 
model with black body. The cutoffpl and the npex 
models showed comparable fit quality and similar sys- 
tematics. As the cutoffpl is the simpler model, we 
will only present these data here. The evolution of 


the spectral parameters is shown in Fig. 6. It is 
clearly seen that neither the photoabsorption Ah nor 
the power law index F seems to vary much over the 
pulse phase. The strong spectral changes indicated 
in the energy resolved pulse profiles are thus mainly 
due to a change in the cutoff energy of the cutoffpl 
model. The cutoff energy is highest in the primary 
peak, making its spectrum distinctly harder than the 
rest of the pulse phase. In the secondary pulse the cut- 
off energy moves to values as low as ~14keV, forcing 
a strong attenuation of the spectrum at higher ener- 
gies. Very interesting is also the behavior of the black 
body component. Its temperature kT is relatively sta- 
ble over the pulse, but its normalization is changing 
by a factor of ~3. The black body is strongest in the 
phase bin between the two peaks and lowest in the 
rise and maximum of the primary peak. This small 
contribution of the black body in the primary peak is 
somewhat compensated by a broadening of the iron 
line in that phasebin. With the energy resolution of 
the PCA it is impossible to tell if the iron line really is 
stronger in that phasebin or if a mixture of blackbody 
and line is fitted. 


V. DISCUSSION & CONCLUSION 

We have presented the first detailed study with IN- 
TEGRAL and RXTE of 4U 1909+07 and have shown 
that the pulse period evolution shows strong indica- 
tion for a random walk like behavior. Such a behavior 
has been seen in other HMXB like Vela X-l and is 
a strong indicator for a wind accreting source with- 
out a persistent accretion disk [18]. An accretion disk 
would provide a more constant transfer of angular mo- 
mentum and thus a long-term spin-up or spin-down 
trend, as seen in other sources like 4U 1907+09 [19]. 
Accretion disks form to fulfill the conservation of an- 
gular momentum of the accreted wind, but in a system 
with a strongly magnetized neutron star the ionized 
matter will couple to the magnetic field lines before 
a stable disk can form. This requires field strengths 
on the order of 10 12 G, not unusual for this kind of 
systems. In these strong magnetic fields, however, cy- 
clotron resonant scattering features (CRSFs) should 
form, as seen in many other sources, e.g., Vela X-l 
[20] or MXB 0656—072 [21], We found no evidence 
for a CRSF in 4U 1909+07. As shown by Schonherr 
et al. [22] this does not rule out a strong magnetic field, 
as CRSFs can be filled up again by photon spawning 
and depend strongly on the geometry of the accretion 
column and on the viewing angle. Further investi- 
gations with high-resolution data will allow a more 
thorough search for a CRSF in 4U 1909+07. A strong 
magnetic field in 4U 1009+07 is also supported by the 
strongly pulsed flux, as the X-ray flux is believed to 
be produced in a very small region above the so-called 
hot-spots on the neutron star surface, where the mag- 
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TABLE I. Fit parameters for the phase averaged spectra with different models 

Model 

Cutoffpl 

Cutoffpl 

Highecut 

Highecut 

NPEX 

NPEX 

parameter 


+bbodv 


+bbody 


+bbodv 

Xred 

1.76 

1.01 

1.07 

0.91 

1.44 

0.93 

Nh [atoms cm* 2 ] 

I5.3t“;« 

4-7 t\i 

4-8lS:S 

a 0+2.0 
0-8 -0.8 

IT 'll 

5-3 

r a 

1-63 1°:£! 

0-96TH 3 

1 Q ’7+0. 03 
1,0 ' -0.08 

1.32 ± 0.10 

0-36:°;°| / -3.08:°;J| 

0.80±g;5| / -2.0 b 

Fe a 

i 7 +oa 

1 -0.5 

o.28±0;if 

n a 1 +0.07 

U.^l_ 0 .o6 

0.2 ±0.2 

n 4i + 0 - 07 
U,4i — 0.06 

0 27 +0 ' 12 
-0.18 

Fe Energy [keV] 

6-4±g;| 


6-39tg;jg 

6.39±g;gt 

a qer+0.04 
o.of _0.03 

6.40 ± 0.05 

Fe Norm 

x io' 3 

! (0.571°; [I ) x 10~ a 

1 (0.78±g;g5) X 10^ 3 

(5±f) x 10~ 4 

(0.80±g;S) X 10- 3 

(0.6±g;f ) x 10~ 3 


a For the NPEX models cvj and ao are given 
15 c*2 was frozen during fitting 



FIG. 6. Parameters of the phase-resolved spectra fitted 
with the cutoffpl model, a ) Pulse profile in the 13.2- 
14.5 keV energy range. The different shaded areas indi- 
cate the phasebins used to extract the spectra, b) Photo 
electric absorption column JVh, c) Power law index F, d) 
Folding energy, e) Black body norm, f) Black body tem- 
perature, g) Iron line width <j. 

netic field lines penetrate the neutron star surface [23] . 
In the accretion column above the hot-spots the den- 
sity and temperatures are highest and thus the most 
X-rays are produced through Comptonization of the 
t hermal photons from the thermal mound on the neu- 


tron star surface [24, 25]. If the magnetic axis is not 
aligned with the rotational axis of the neutron stars, 
the hot spots and the accretion columns move in and 
out of sight of the observer, resulting in a pulsed X-ray 
flux. 

The two distinct peaks in the pulse profile at low 
energies hint at the fact that accretion happens onto 
both magnetic poles of the neutron star, but under 
different physical conditions. The different conditions 
can explain the different spectra in the two peaks. The 
secondary peak shows a distinctly lower cutoff energy 
which corresponds roughly to the temperature of the 
electron gas in the accretion column. A possible ex- 
planation could be, e.g., different sizes of the hot-spots 
due to misalignment of the magnetic center and the 
gravitational center of the neutron star. This would 
lead to different sizes in the accretion column and thus 
to different densities and temperatures. 

The exact geometry of the accretion column de- 
pends strongly on the way the matter couples to the 
magnetic field and is highly uncertain [26. 27]. Besides 
the simple filled funnel other possible configurations 
include a hollow or a partly hollow funnel. With an 
average luminosity of ~ 2.8 x 10 36 ergs -1 the system 
can host accretion rates that are large enough for a 
shock to form in the accretion column, due to the 
radiation pressure from the material close to the neu- 
tron star surface. This shock prevents escape of the 
Comptonized photons parallel to the magnetic field 
lines, so that they can only escape on the sides of the 
accretion column in a so-called "fan beam” [27]. If the 
accretion rate is lower, no shock will form and the ra- 
diation can escape along the accretion column in the 
"pencil beam” . Which of these two geometries exists 
in a given object is hard to tell from timing analysis 
alone, as no independent measurements of the align- 
ment of the magnetic axis is possible. But phase- 
resolved spectroscopy can provide some insight into 
the configuration. For 4U 1909+07 we have seen that 
the blackbody component of the spectrum is phase 
shifted to the two peaks in the pulse profiles. It is most 
likely that the blackbody originates from the ther- 
mal photons of the thermal mound at the hot spot. 
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Through the rotation of the neutron star, the hot 
spot is only visible at certain pulse phases. Assum- 
ing a hollow funnel accretion geometry, with a large 
enough accretion rate to form a shock, we can argue 
that we view the thermal mound only at the phases 
were we look along the accretion column, while t he 
most hard X-ray flux is only visible when we see the 
sides of the accretion column, due to the “fan beam” 
geometry. 

Caution has to be taken, however, with these mod- 
els as the behavior of plasma in magnetic fields on the 
order of B as 10 12 G is still only little understood and 
all models are based on gross simplifications. Addi- 
tionally gravitational light bending must be taken into 
account when analyzing pulse profiles and their origin 
[28] . Detailed analysis of this kind is ongoing and will 
be presented in a separate publication. 
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